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Abstract. Critical states are studied by a real space RG in the problem with strong diagonal 
disorder and long range power law hopping. The RG flow of the distribution of coupling 
parameters is characterized by a family of non-trivial fix points. We consider the RG flow 
of the distribution of participation ratios of eigenstates. Scaling of participation ratios is 
sensitive to the nature of the RG fix point. For some fix points, scaling of participation ratios 
is characterized by a distribution of exponents, rather than by a single exponent. 

The RG method can be generalized to treat certain fermionic Hamiltonians with disorder 
and long range hopping. We derive the RG for a model of interacting two-level systems. 
Besides couplings, in this problem the RG includes the density of states. The density of 
states is renormalized so that it develops a singularity near zero energy. 
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Conventionally, localization is studied for Hamiltonians virith matrix elements rapidly 
decreasing as a function of the distance. For instance, the canonical Anderson model is 
defined on a lattice so that disorder is diagonal and hopping occurs between neighbor- 
ing sites. In this model the localization transition is reached by reducing the amplitude 
of hopping V below the threshold Vc approximately given by the disorder bandwidth 
W divided by the number of nearest neighbors. 

The situation is completely different for the problem with long range hopping|]lj. 
In a d— dimensional space, for hopping amplitude falling as some power of the distance, 
V{r — r') ^ |r — r'|~", all states are extended for a < d for both strong and weak 
disorder. At the same time, for a > d and strong disorder the states are localized. 

Thus, one can study localization transition by varying the exponent a at fixed 
disorder strength ^, ^, ||] . This transition is more tractable than the conventional 
one because of the availability of a small parameter given by the ratio of hopping to 
disorder strength At the transition line a = d, a real space renormalization 

group (RG) can be constructed for the distribution of couplings. The RG flow obeys 
certain conservation laws, and, depending on the details of the Hamiltonian, the flow 
may or may not have a nontrivial flx point. 

Recently, this problem reappeared in the context of quantum chaos in Kepler 
billiards The RG approach was used in [|| to treat scaling of wave functions at 
the transition. In Sec. ^ below the RG scheme is applied to study the distribution of 
participation ratios. We find that, depending on the type of the fix point, participa- 
tion ratios are either characterized by a single scaling exponent, or by a distribution 
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of scaling exponents. 

One can extend the RG method to certain many-body systems with disorder 
and long range interactions. One interesting example is the quantum problem of two 
level systems with long range coupling |§, 0, § . In Sec. |^ we develop an RG for a c? = 3 
problem with couplings falling as a cube of the distance. From the RG flow derived 
for the distribution of the wavefunctions and for the density of states we find that: 

(i) Wavefunctions of all excitations are delocalized; 

(ii) Under RG of the density of states a singularity develops near zero energy due to 
states with small energy condensing into a peak. We discuss the relation of this new 
state of interacting two level systems to the Efros-Shklovskii Coulomb gap state M. 



1 Renormalization group 

We consider an ensemble of Hamiltonians on a d— dimensional lattice: 

E-ipr = Er^pr + V Vr^r'V'r' , whcrC Vr,r' = i"'' "■'"^ ■ (1) 

Here are random numbers, and are random n— component vectors. The distri- 
bution of Ej. is uniform with density v. dN/dE — v for \E\ < l/2z/; for \E\ > l/2i/. 
We assume that the distribution of dr, dP = f{d)d'^a, is such that all moments of 
/(a) exist. Otherwise, the exact form of /(a) is unimportant. 

The coupling strength is characterized by dimensionless parameter A = {a1)nv, 
where n is the density of lattice sites. In the weak coupling regime, A ^ 1, a real 
space RG flow of the Hamiltonians of the form (|l|) can be derived. This is done by 
considering resonance pairs and arguing that they correctly account for the flow of the 
coupling distribution to the leading logarithmic order (see !§,§]). 

The RG flow preserves the form of the hopping amplitudes Vi-^r' ui (^, and modifies 
the distribution /(a). The distribution of energies E^ is renormalized near the tails 
but remains constant in the middle part. In the middle of the spectrum, |£'r| <C , 
one can treat Ej. as quasirandom numbers with fixed uniform distribution. 

To introduce some objects appearing in the RG equation, let us explicitly diago- 
nalize the problem of a resonance between two sites ri and V2'. 

ElPi = Ei^i + ^ri,r.V'2 , SV^2 = -^2^2 + 14,,riV'l • (2) 

Two eigenstates tj}^ and are defined by 

V'^ = cos^ Vi + sin^ V2, -0" = ~sin0 -|- COS0 V2 , (3) 

with 

tan 61= Vr^ + l-r , ^ ^ Ei - E2 

2Vi-j^ ^1-2 

The energies of the states (||) are 
1 



E± = -[E,+E2± ^liE^ - E2Y + ) (5) 
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The transformation rule for the vectors 01^2 foUows from the relation 

aiipi + a2'4'2 = a+if}'^ + a-ip~ . (6) 
One obtains 

a+ = cos 9 ai + sin 9 0,2 , a_ = — sin 9 di + cos 9 (12 (7) 

The role of transformed d± is that they determine hopping matrix elements Vi-,r± = 
Or • d±/\v — v±\'^ for all remote sites |r — r±| » |ri — r2|, where r± — (ri + r2)/2. 

Subsequently diagonalizing pair resonances at all distances, from small to large, 
one derives an RG equation for the distribution of a^s. It has a form of a first 
order differential equation with respect to the RG time i = In [L) with an integral 
operator on the right hand side: 
d 



Q-fia) = 2Adi^n J dTd'^aid'^a2 |ai -021 /(ai)/(a2) 

(5 (a — a+) + (5 (a — a_) — (5 (a — ai) — 5 {a — a2) . (8) 

Here ^d, the area of a d— dimensional sphere, arises from dV = A^r'^^^dr. 

The flow of the distribution f {a,t) obtained by solving (||) has a number of prop- 
erties similar to that of Boltzmann kinetic equation. The reason is that in considering 
only pair resonances and ignoring correlations between resonances at different scales 
one makes an approximation of the same nature as the pair collision approximation 
in Boltzmann theory of weakly nonideal gases. 

Two properties of the flow M) are of interest. First, all second moments 



Gap = {aaap) = J f{a)aaapd^a (9) 

are conserved. Second, one can define entropy 

H[f] = - J ln(|a|/(a))/(a)rf"a (10) 

and, for the flow (||) modified by |ai • 02! — > |ai||a2|, prove an H-theorem, dH/dt > 0. 

These results can be used to fully characterize the fix points of the flow (||) . By 
maximizing the entropy (^o|) for given moments Gap one obtains the family 

fcia) ^ {A/\a\)exp{~aaGapap) . (11) 

The normalization constant A depends on the matrix G. 

It is interesting that the distributions ([ll] ) are normalizable only when the number 
of components n is greater than one (and matrix Gap is nondegenerate) . Because of 
that, the flow (^ has very different properties for n = 1 and n > 1. 

In the n > 1 case, for any initial conditions f{a)t=o, the distribution /(a) evolves 
to one of stationary distributions (pT|). On the other hand, for n = 1 possible fix 
point distributions allowed by H-theorem are non-normalizable. Therefore, since 
/ f{a)t=oda = 1, the flow (H) does not have a fix point. 

These two qualitatively different types of flow are displayed in Figures |^,^. The 
simulation was performed by using TV — 10* states characterized by random energies 
and coupling vectors. The energy distribution was uniform with v = I and did not 
evolve. The distribution of a evolved to a fix point of the form ( pi] ) for n — 2 and 
diverged for n ~ I. 
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Fig. 1 The flow (g) for n — 2. Initial distribution is uniform in tfie dislc + a2 < r^, and 
hence f(a) is circular symmetric. At large times the distribution approaches a fix point of 
the form ( |ll| ) Parameters used in this simulation are: v = 1, ro = 0.6, T = 4. Dimensionless 
coupling A = 0.2 for these parameters. 
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Fig. 2 The flow (^ for n = 1. Initial distribution is uniform in the interval \a\ < ro- At 
large times the distribution looks more and more like {A/\a\)e~^°' with time-dependent A 
and g. Note the divergence developing near a = and in the tails. Parameters used in this 
simulation are: f = 1, ro — 0.5, T = 4. Dimensionless coupling A = 0.08. 
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2 Scaling of wavefunctions 

The RG hierarchy of states constructed as resonances on different scales implies fractal 
structure the wavefunctions. To understand scaling of eigenstates, one can extend RG 
by including in it the distribution of wavefunction amplitudes. 

Below we consider the distribution of participation ratios of the states. For a 
normalized state V'r'' its participation ratio p4 is defined as p4 = ^ |V'r*''|* In the 

r 

localization problem the participation ratios are used as a measure of the degree of 
localization. The scaling exponent fi of participation ratios defined as p ^ L~'^, where 
L is the system size, characterizes possible universality classes. In the localized regime 
H = 0, in the delocalized regime n = d, and at the critical point < jj, < d. 

It is straightforward to put the distribution of participation ratios in the RG 
scheme. For that one has to consider the flow of the distribution f(a,p). Follow- 
ing the discussion of Sec. |l|, one finds the participation ratios p± of the eigenstates of 
the resonance pair (H) in terms of pi and P2- Since the states 1 and 2 do not overlap, 
one simply obtains: 

p+ = cos^ d pi + satl^ 6 p2, P- — srn^ Pi + coB^ 6 p2. (12) 

This change of the participation ratios must be considered together with the transfor- 
mation of the parameters oi, 02 given by Eq. (|^). The resulting RG equation reads: 

d f 

— f{a,p) = 2Adm / drd" aid"- a2dpidp2 |ai-a2| / (ai,pi) / (02,^2) 

X [S (a - a+) 5 [p - p+) + 5 {a — a-) S {p - p_) 

-d{a-ai)S{p-pi)-dia-a2)S{p-p2)] . (13) 

Generally, one cannot factor the distribution f{a,p), because the flow ( p^ ) generates 
nontrivial correlations between a and p. 



It is of interest to consider the flow (13) in the scaling limit, when the distribution 
of a given by J f{a,p)dp evolves to the fix points discussed in Sec. |l|. We studied 
numerically evolution of the participation ratios distribution J f{a,p)d^a. 

Because of the expected power law scaling p ^ = exp(— /xi) it is natural to 
consider the distribution of In p. The results of the simulation for n — 2 and n — \ are 
shown in Figures Parameters used are the same as in Figures |^J|. 

The results for n — 2 and n — 1 are quite different. For n = 2, the distribution 
of \np forms a peak of slowly varying width. The center of the distribution moves 
linearly with the RG time which corresponds to a single scaling exponent /i given 
by the peak velocity. In contrast, for n = 1, the left edge of the distribution is fixed at 
hip = 0. After rescaling, the distributions at different t collapse on one another. This 
means that participation ratios are characterized by a continuous spectrum of scaling 
exponents /i, rather than by a single exponent. 

The values of /i are not universal because of a many-parameter family of fix points 
parameterized by Gaf3 as discussed in Sec. |^. Even for spherically symmetric distri- 
butions with Gap = Gdaf}, the distribution depends on the dimensionless parameter 
A — Gnv. The dependence on A however is quite simple. Numerically, we find that 
doubling A doubles velocity of the center of the distribution of hip, without changing 
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of its overall form. This implies that scaling exponents /i can be written as 7A, where 
7 ~ 1 and has a universal distribution. 

This estimate is in agreement with the result for the mean participation ratio [§, 
where it was found that (p) scales as L~'^^ with 7 ~ 1. 



3 RG for fermionic Hamiltonians 

In some cases, the RG approach outlined above can be extended to many-body prob- 
lems Here we consider a Hamiltonian for interacting two-level systems in d = 3: 

w = E + ^ E '^i^'^^f ^"'^^ (14) 

Here af are pseudospin Pauli matrices (a, b — x,y, z), and r.y = jr.; — rj \ are distances 
between two-level systems randomly distributed in space with concentration n. Both 
ei and V^"'' are random uncorrelated numbers. 

The problem ( |l4|) describes several physical situations, such as the dynamics of 
excitations in two-level systems in glasses interacting via elastic strain or dipole 
excitations in the Efros-Shklovskii model of a disordered Mott insulator [Q. In the 
latter problem, the conventional approach is to consider only crf(T| couplings in (^4|). 
In this case excitations are strictly localized, and the only effect of interaction is the 
appearance of a logarithmic gap in the density of states at low energies. 

Here we focus on the opposite limit, and consider the effect of x — x, y — y, and 
X — y couplings, ignoring other couplings. We assume that the interaction is weak, 
A = ^ 1, and develop an RG scheme using this weak coupling. The 
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participation ratio, In(P^) 



Fig. 4 The flow of participation ratios for n = 1 corresponding to the parameters used in 
Figure |[ The RG times = 2mT, 1 < m < 10. 



XY approximation is frequently used for studying spin waves in spin systems. Al- 
though it does not seem to be controlled by any small parameter, one should mention 
that the terms removed have zero matrix elements between different states with close 
energy. Therefore, these terms cannot produce resonances (which are responsible for 
delocalization), so they are unlikely to be important as long as delocalization of states 
is concerned. 

Fermionic model is obtained by replacing Pauli matrices in ( [T^ by Fermi operators: 
—>■ a + a'^ , ct^ — > i{a — a~^), a' — > a^a — aa^. Our main reason to employ this spin- 
pseudofermion transformation is that it gives an exact result for a pair of interacting 
spins, whereas only interacting pairs contribute to the leading logarithmic order in the 
RG. Thus we obtain 

H = '^ei{alai - aiol) + ]^'^^r^j^{uijajaj + Vijaittj + H.c.) , (15) 

i i=£j 

where tty = + V^f + i(l^,7 - V^f ), = V^f - Vf^ - i{V^^ + V;f ). Distributions 
of real and complex Uij , Vij are assumed to be uncorrelated. We denote them by vU) 
and /(w, f), respectively: dP = v{e)de and dP = f{u, v)d'^ud'^v. The Hamiltonian ( [l5| ) 
is bilinear in a, o"*" and can be treated perturbatively by an RG in the weak coupling 
regime {\u\)vn<til, (|w|)w<Cl. 

The RG scheme is based on resonance pairs of pseudofermions: 

Ti = ei{d^ai — aia^) + £2(02^02 ^ 02^2^) + r^2{uaia2 + waia2 + H.c ) (16) 
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The Hamiltonian ( ^q ) can be diagonalized by a two step Bogoliubov transformation: 
(1) (0,1,0,2) = Ra{ai,a2); (2) ia'i,a'2) = Rpia^, 02), where 

_ ( cos a zisinaX Ji ^ ( '^^^(^ sin /?\ 
" \^ — sin a zi cos a J ' ^ \^ — sin /? Z2 cos /? / ' 

— 3 —3 

tan2a = ^^ — , zi = —, tan2/3 = — , Z2 = ri ■ 18) 

£1-62 \u\ £1+62 |w| 

The resuhing diagonal Hamiltonian is 7-^ = e'i{a'^a'i — a'la'^) + £2(^12^02 ^ 02*^2^)1 
where 

ie\ - e'2f = (£1 - £2)^ + \ur^^\\ ie\ + ^'2? = (£1 + £2)' + \vr^^?. (19) 

For another two-level system (say, described by b and 5+) the interaction with ai and 
02 is given by 

'Hmt ~ r^^ [uio^b + viaib + H.c.) + {u2a2b + V2a2b + H.c.) (20) 

In accordance with the picture of resonance pairs we can consider only the case 
''ib^'^12, T2b^fi2- Therefore, rib is almost equal to r2b and hence the two Bogoliubov 
transformations act on U1.2 and vi,2 as 

({ti,U2) = Ra{ui,U2), {vi,V2) = i?* (t)l , ^2 ) ; 

(Ui,W2) = i?^(Ml,'52), (1)^,^2) = i?/3(wi,U2) , (21) 

where * stands for complex conjugate. Similarly to how it was done in Sees |^J^, these 
expressions can be used to derive an RG flow for distributions: 

|g(0 = 2^ne''m ei) + ^{i - Q - - a) - ^(e - , (22) 

Here we introduced a condensed notation = (ui,Vi,ti), so g(^) is a distribution of 
M, w, and £, the 5— function — — 5{u — ui)5{v — vi)5{e — £i), etc. The brackets 
(■■■)« t) 5i ?2 st^'^cl for averaging J ...f{u,v)(Pu(Pvg{^i)d^^ig{S,2)d^S,2- (Here the inte- 
gration volume elements are 2- and 5-dimensional because u and v are complex.) The 
expression (...)„ „ in ( p2| ) is of order of e"''', so the RHS of ( |2^ ) has no explicit 
dependence on t. 

The RG flow (p2|) restricted to the distribution of couplings u, v is analogous to 
that discuss in SecTyJ. One can prove an H-theorem, verify that second moments are 
conserved, and derive stationary fix-point distributions: 

fw,z{u, v) = TT^'^{w'^ - |zp) exp(-i(7(|up + |up) - {zuv + c.c.)) (23) 

where w is real and z is complex, satisfying \z\ < w for normalizability of fw,z{u, v). 

A new aspect of the pseudofermion problem is a nontrivial RG flow of the density 
of states. To study how the distribution of energies e evolves, we consider a particular 
class of distributions: 

g{u, V, e) = iy{e)fw,z{u, v) (24) 
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with /iu.z(u,w) being the fix-point distribution of the form (p3|). (Such an ansatz is 
consistent with Eq.(p2|).) Substituting this into ( p^ ) one obtains 

^Ke) = 2nne''{d{e - e[) + S{e ~ e',) 6{e e,) - 6{e - 62))„,„,,„,, (25) 

We studied the flow of t^(e) numerically (see Figure ||). There are two distinct features 
in the flow of v{e). Initially, at small t, the density of states decreases at all e due to 
energy repulsion. At large i, it increases in the vicinity of e = forming a singularity. 
It is clear from Figure |^ that the first effect is transitional and quickly saturates, 
whereas the second one is characteristic for the scaling limit. 

It is interesting to compare this result with the Efros-Shklovskii theory of Coulomb 
gap which predicts a logarithmic suppression of the density of states near e = for the 
problem ( p^ with z — z couplings. In the above analysis we ignored these couplings 
because they do not lead to delocalization of excitations. The effect of 2; — 2 couplings 
on the density of states is of the same order as the effects we consider. However, it 
has an opposite sign. 

Therefore, there are two different phases possible in the problem (|l^ depending 
on relative strength of different couplings. If the z — z coupling dominates, the system 
will have a soft gap in the density of states. On the other hand, if the x — x, y — y, 
and X ~ y couplings dominate, there will be a peak in the density of states at e = 0. 
Our analysis indicates that the peak is described by a power law singularity with the 
exponent of order of A. 
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In both phases, however, excitations are delocahzed. According to Sees. |l|,y|, the 
states in this problem are critical. Delocahzed states will give rise to thermal conduc- 
tivity and, in the Coulomb gap state, will make variable range hopping conductivity 
possible even in the absence of electron-phonon coupling. 



4 Summary 

The RG approach predicts nontrivial fix points of Hamiltonians with long range hop- 
ping. Wavefunctions in such systems are critical, and the participation ratios are 
characterized by scaling exponents which depend on the nature of the fix point. In 
some cases, for a particular Hamiltonian different states can have different exponents, 
and so the system is characterized by a distribution of exponents, rather than by a 
single exponent. (The possibility of the absence of self-averaging of the participation 
ratio exponent was considered recently pl|.) 

The RG method can be applied to the problem of interacting two-level systems. 
In doing this it is important to extend the RG by including renormalization of the 
density of states. Depending on the choice of couplings, two different phases are 
possible, both characterized by a singularity in the density of states at e = 0. In one 
phase the density of states diverges at e = 0, whereas in the other phase it vanishes. It 
is possible that the first behavior is relevant for the problem of low energy excitations 
in glasses, where an increase in the density of states at low energies due to interactions 
between two-level systems has been conjecturedpOf. 

I benefited from discussions with B. L. Altshuler, V. E. Kravtsov, A. D. Mirlin, and B. I. Shklovskii. 
This research is supported in part by the MRSEC Program of NSF under award 6743000 IRG. 
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